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HIGHLIGHTS 


•  Performance  degradation  of  lithium-ion  battery  over  the  cell  lifetime  is  quantified. 

•  State  estimators  with  the  different  time  scales  are  developed  for  SOC  and  SOH  identification. 

•  Capacity  fade  and  power  fade  are  accurately  characterized. 

•  SOC  estimator  is  accurate  and  robust  over  the  life  span  of  the  battery  cell. 
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A  combined  SOC  (State  Of  Charge)  and  SOH  (State  Of  Health)  estimation  method  over  the  lifespan  of  a 
lithium-ion  battery  is  proposed.  First,  the  SOC  dependency  of  the  nominal  parameters  of  a  first-order  RC 
(resistor-capacitor)  model  is  determined,  and  the  performance  degradation  of  the  nominal  model  over 
the  battery  lifetime  is  quantified.  Second,  two  Extended  Kalman  Filters  with  different  time  scales  are 
used  for  combined  SOC/SOH  monitoring:  the  SOC  is  estimated  in  real-time,  and  the  SOH  (the  capacity 
and  internal  ohmic  resistance)  is  updated  offline.  The  time  scale  of  the  SOH  estimator  is  determined 
based  on  model  accuracy  deterioration.  The  SOC  and  SOH  estimation  results  are  demonstrated  by  using 
large  amounts  of  testing  data  over  the  battery  lifetime. 
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1.  Introduction 

Lithium-ion  batteries  have  been  widely  used  in  modern  elec¬ 
trified  vehicles.  The  reliable,  efficient,  and  safe  operation  of  lithium- 
ion  batteries  requires  monitoring,  control  and  management.  For 
battery  management  systems,  a  core  function  is  to  provide  accurate 
estimates  of  State  of  Charge  (SOC)  and  State  of  Health  (SOH)  of 
batteries,  which  is  challenging  due  to  the  lack  of  sensors  for  elec¬ 
trochemical  phenomena  inside  the  cells. 

Many  methods  were  proposed  to  estimate  the  battery  SOC,  each 
with  its  own  advantages  and  disadvantages,  as  summarized  in 
Table  1.  The  Coulomb  counting  method  and  open  circuit  voltage 
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method  are  widely  used  in  battery  management  systems  of  elec¬ 
trified  vehicles.  They  are  easy  to  use  and  fast  in  computation,  but 
the  former  highly  relies  on  the  performance  of  current  sensor,  and 
the  latter  is  not  effective  for  batteries  with  flat  open-circuit-voltage 
curve.  Another  disadvantage  for  Coulomb  counting  is  that  this 
method  is  open-loop  estimation  and  may  have  large  accumulated 
error  due  to  uncertainties  or  disturbances  [1—5],  Moreover,  it  re¬ 
quires  accurate  initial  SOC  value.  Many  artificial  intelligence-based 
methods  have  been  applied  to  establish  black-box  SOC  estimation 
models,  such  as  neural  network  [6],  fuzzy  logic  [7],  and  support 
vector  regression  (SVR)  models  [8],  The  Kalman  filter  (KF)  and 
sliding  mode  observer  have  also  been  used  to  predict  the  battery 
SOC.  These  approaches  are  model-based,  closed-loop,  and  thus  can 
use  output  feedback  to  keep  better  robustness  than  non-feedback 
methods.  In  Refs.  [9—11],  the  extended  Kalman  filter  (EI<F) 
concept,  based  on  nonlinear  state-space  models,  was  used  to 
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Table  1 

Advantages  and  disadvantages  of  existing  SOC  and  SOH  estimation  methods. 


State  of  charge  (SOC) 

State  of  health  (SOH) 

Method 

Advantage 

Disadvantage 

Method 

Advantage 

Disadvantage 

Coulomb  counting 

Simple 

Open-loop,  sensitive  to  the 

Durability 

Durability  mechanism 

Comprehensive 

Complex,  need 

[1-4] 

current  sensor  precision,  and 
uncertain  to  initial  SOC 

model-based 

open-loop 

[23,24] 

understanding 

accurate  input 
parameters 

Open  circuit 
voltage 
method  [5] 

Simple 

Open-loop,  sensitive  to  the 
voltage  sensor  precision, 
unsuitable  for  cells  with  flat 
OCV— SOC  curves 

method 

Durability  external 
characteristic  [25—28] 

Simple  and  easy  to 
predict  capacity 
fade  and  internal 
resistance  increment 

Based  on  a  large 
number  of 
experiments 

Neural  network  [6] 

Generic,  good  nonlinearity 
mapping  approximation 

Sensitive  to  the  amount  and 
quality  of  training  data 

Battery 

model-based 

parameter 

DC  resistance  [21] 

Simple 

Not  accuracy, 
sensitive  to 
disturbances 

Fuzzy  logic  [7] 

Generic,  good  nonlinearity 
mapping  approximation 

Sensitive  to  the  amount  and 
quality  of  training  data 

identification 

closed-loop 

AC  impedance  [22] 

Accuracy 

Complex 

Support  vector 

Generic,  good  nonlinearity 

Sensitive  to  the  amount  and 

method 

Extend  Kalman  filter 

Quite  easy  to 

Sensitive  to  modeling 

machine  [8] 

mapping 

quality  of  training  data 

[11,29] 

implement,  accurate 

accuracy 

Kalman  filter 

[9-18] 

Closed-loop,  online, 
accuracy 

More  computationally 
expensive  than  non-feedback 
methods,  and  highly  depend 
on  the  model  accuracy. 

Fuzzy  logic  [30] 

Sample  entropy  [31—33] 

Accuracy  simple, 
accurate 

Simple 

Slow  convergence 

Need  large  amount 
of  data 

Sliding  mode 

Closed-loop,  online,  and 

More  computationally 

Discharge  voltage  [30] 

Easy 

Not  accurate 

observer 

[19,20] 

accurate 

expensive  than  non-feedback 
methods,  and  highly  depend 
on  the  model  accuracy. 

Adaptive  control  system 
[31] 

Online 

Sensitive  to  modeling 
accuracy 

estimate  the  SOC  of  a  Li-polymer  battery.  Several  other  variants  of 
Kalman  filter,  e.g.  sigma-point  KF  [12,13],  adaptive  KF  [  14—16],  Dual 
KF  [17]  and  derivative  KF  [18]  have  also  been  used  for  battery  SOC 
estimation.  The  sliding-mode  observer  technique  has  also  been 
used  to  monitor  battery  SOC  trajectories  [19,20], 

State  of  Health  (SOH)  is  a  metric  to  evaluate  the  aging  level  of 
batteries,  which  often  includes  capacity  fade  and/or  power  fade. 
The  commonly  used  indicators  include  battery  capacity  [11],  DC 
resistance  [21],  and  AC  impedance  [22],  The  SOH  estimation 
methods  mainly  include  durability  model-based  open-loop 
methods  and  battery  model-based  closed-loop  method  [5].  The 
former  methods  directly  predict  the  changes  in  capacity  fade  and 
internal  resistance.  The  durability  models  describe  the  increase  of 
SEI  film  resistance  and  battery  terminal  voltage  [23,24],  Based  on 
durability  characteristics,  a  storage  life  model  for  lithium  cobalt 
oxides  batteries  was  given  in  Ref.  [25].  Bloom  et  al.  [26]  obtained 
the  relationship  between  the  battery  performance  degradation  and 
ambient  temperatures  and  cycle  time.  Matsushima  [27]  also  found 
that  capacity  loss  exhibits  a  square  root  relationship  with  time.  Li 
et  al.  [28]  developed  an  extended  Arrhenius  model.  The  battery 
model-based  closed-loop  methods  use  least-squares  methods, 
Kalman  filtering  [29]  and  other  adaptive  algorithms  (such  as  fuzzy 
logic  [30]),  to  identify  the  battery  capacity  and  internal  resistance 
according  to  the  operating  data.  Sample  Entropy  was  also  used  to 
estimate  the  battery  SOH  in  Refs.  [31—33].  The  advantages  and 
disadvantages  of  these  SOH  methods  are  summarized  in  Table  1. 

Most  of  the  above  mentioned  battery-state-estimation  methods 
were  developed  for  either  SOC  or  SOH  estimation,  and  not  both.  The 
intimate  coupling  feature  between  SOC  and  SOH  was  overlooked. 
The  accuracy  of  SOC  estimation  is  heavily  influenced  by  battery 
degradation.  As  batteries  degrade,  SOC-only  estimation  algorithms 
may  lead  to  large  errors.  The  inaccurate  SOC  estimations  in  turn 
may  mislead  the  battery  SOH  calibration.  Therefore,  simultaneous 
estimation  of  SOC  and  SOH  is  quite  beneficial.  Compared  to  the 
battery  SOC  variation,  battery  SOH  typically  change  much  more 
slowly,  necessitating  multi-timescale  state  estimators.  In  order  to 
determine  the  appropriate  time  scale  for  the  SOH  estimator,  it  is 
critical  to  examine  the  performance  degradation  of  the  battery 
model  in  the  context  of  battery  aging.  The  multi-scale  EKFs  are  used 


to  estimate  SOC  and  SOH,  and  the  capacity  estimation  is  periodi¬ 
cally  introduced  in  SOC  update  equation  [34,35].  It  is  more 
computational  efficient  than  a  joint  estimation  [34],  However,  the 
determination  of  the  two  time  scales  is  heavily  dependent  on  the 
tuition  and  calibration. 

This  paper  discusses  a  model-based  combined  SOC/SOH  esti¬ 
mation  method  over  the  lifespan  of  LiNMC  batteries.  First,  the  SOC 
dependence  of  the  nominal  parameters  of  a  first-order  RC  (resis¬ 
tor-capacitor)  model  is  determined,  and  the  performance  degra¬ 
dation  of  the  nominal  model  is  quantified  over  the  battery  lifetime. 
Second,  two  EKFs  with  different  time  scales  are  applied  to  imple¬ 
ment  the  combined  SOC/SOH  monitoring:  one  observer  is  for  real¬ 
time  SOC  estimation;  the  other  for  offline  SOH  (capacity  and  in¬ 
ternal  resistance)  update.  The  time  scale  of  the  SOH  estimator  is 
determined  based  on  the  quantified  model  accuracy  degradation. 
The  SOC  and  SOH  estimation  results  are  demonstrated  by  using 
large  amounts  of  testing  data  over  the  battery  lifetime. 

The  remainder  of  the  paper  is  arranged  as  follows:  Section  2 
introduces  the  battery  model  structure  and  battery  tests;  the 
degradation  of  a  nominal  battery  model  is  described  in  Section  3; 
the  combined  SOC/SOH  estimation  approach  and  associated  results 
are  discussed  in  Section  4;  Section  5  concludes  this  paper. 


2.  Battery  modeling 

2.1.  Equivalent  circuit  model 


Hu  et  al.  [36]  compared  12  commonly  used  equivalent  circuit 
models  and  concluded  that  the  first-order  RC  model  is  the  best 
choice  considering  model  complexity,  accuracy,  and  robustness. 
Fig.  1  shows  the  model  structure  considered.  The  battery  capacity 
Ccap  is  used  to  quantify  SOC  level  by  Eq.  (1) 


SOC=  - 


V '  A>att 

3600  Ccap 


(1) 


The  Coulombic  efficiency  r/  is  simplified  as  the  constant  value, 
1.0  during  the  discharge  and  0.98  in  charging. 
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SOC  Description  Voltage-Current  Characteristics 


Fig.  1.  Battery  model  structure:  open  circuit  voltage  Voc,  series  resistance  Ro,  diffusion 
resistance  Rp,  and  diffusion  capacitance  Cp;  current  /batt.  terminal  voltage  Vt,  and  bat¬ 
tery  capacity  Cap- 


The  capacity  calibration  is  carried  out  by  Constant  Current 
Constant  Voltage  (CCCV)  reference  test.  Controlled  current  source 
/batt  representing  the  current  flow  in  the  cell,  acts  as  the  input 
variable.  Controlled  voltage  source  Voc  stands  for  the  cell  Open 
Circuit  Voltage  (OCV).  Terminal  Voltage  Vt  is  the  output  variable. 
The  ohmic  resistance  Ro  describes  the  internal  ohmic  resistance. 
The  mathematical  equations  are  shown  below: 

fup  =  -3v  +  /^tt 

<  cp^p  CP  (2) 

1  Vt  =  Voc  —  Up  —  /batt^O 

where  up  is  the  voltage  across  the  parallel  RC  network. 

2.2.  Battery  testing  systems  and  schedule 

The  test  data  used  for  this  study  were  acquired  through  the 
battery  testing  facility  in  Refs.  [33,36],  The  data  acquisition  fre¬ 
quency  is  10  Hz.  The  batteries  used  for  this  test  are  lithium  nick¬ 
el— manganese— cobalt  (LiNMC)  UR14650P  cells  from  Sanyo  with  a 
graphite  anode.  Its  main  specifications  are  listed  in  Table  2.  Eight 
cells  were  tested  with  the  compound  load  profiles  for  comparison 


reached.  More  detailed  discussion  of  the  battery  testing  systems 
and  schedules  can  be  found  in  Ref.  [36], 


3.  Battery  model  parameter  identification 

3.1.  Identification  method 


Fig.  3  shows  the  data  collected  in  the  hybrid  pulse  test  at  22  °C, 
and  Fig.  4  is  a  magnified  view  of  one  discharge  pulse  at  around  90% 
SOC.  The  detailed  procedure  is  illustrated  as  follows: 


(1)  The  Voc  value  is  determined  from  the  steady-state  voltage,  as 
shown  in  Fig.  4(b). 

(2)  The  ohmic  resistance  R0  is  proportional  to  the  instantaneous 
voltage  drop  dV  following  a  discharge  pulse. 


(3) 


(3)  The  polarization  resistance  and  capacitance  (Rp,  Cp)  are 
related  to  the  voltage  amplitude  change  following  the 
instantaneous  potential  variation  described  above,  as  shown 
in  Fig.  4(b).  The  Recursive  Least  Squares  (RLS)  algorithm  is 
used  for  (Rp,  Cp)  extraction.  The  linear  identifiable  form  is 
derived  by 


-  V(k)  -  Voc  (k)  +  Voc(/<  -  1)  -  V(k  -  1)  =  ai/(k)  +  a2I(k  -  1) 

-RpCpR0  —  (Rp  +  Ro)  At 

-  “1= - At - 

RPCpRo 

(4) 

where  a ^  and  a2  are  two  parameters  to  be  identified  by  RLS,  given 
Voc  and  Ro  from  the  above  two  steps.  The  associated  RLS  parame¬ 
terization  is  described  by 


gw  =  y-wQ  , 

1  +  </>T(k)P(k-  \)<p(k) 

e(k)  =  0{k  -  1)  +  G(k)  [V(l<)  -  Voc (k)  +  Voc (k  -  1)  -  V(k  -  1)  -  </>T(k)9(k  -  1)] , 
-  P(k)  =  P(k  -  1)  -  G(k)vT(k)P(k  -  1), 


and  validation.  As  shown  in  Fig.  2,  the  experiment  procedure  starts 
with  a  series  of  characterization  tests  at  three  different  tempera¬ 
tures  (in  the  order  of  10  °C,  35  °C  and  22  °C).  Each  characterization 
test  consists  of  a  static  capacity  test  using  current  rate  0.5C,  a  hybrid 
pulse  test,  a  DC  resistance  test,  a  dynamic  stress  test  (DST),  and  a 
Federal  Urban  Driving  Schedule  (FUDS)  test.  After  these  tests,  aging 
cycles  are  conducted  at  22  °C.  In  each  aging  cycle,  the  cells  are 
charged  and  discharged  at  a  constant  rate  until  the  cut-off  voltage  is 


Table  2 

Main  specifications  of  the  LiNMC  cell. 


Nominal 

Nominal 

Upper  cut-off 

Lower  cut-off 

capacity  (Ah) 

voltage  (V) 

voltage  (V) 

voltage  (V) 

0.94 

3.70 

4.20 

2.50 

where 

(cp(k)  =  m,i(k--i)}T,  (6) 

[9(k)  =  [a-l,a2]r. 

The  initial  values  of  the  parameter  estimate  0(0)  and  its  error 
covariance  matrix  P( 0)  need  to  be  firstly  prescribed.  After  cali¬ 
brating  9,  Rp  and  Cp  can  be  easily  determined  by  Eq.  (4).  The  stan¬ 
dard  RLS  algorithm  and  some  variants  were  often  used  for  battery 
modeling  [37,38], 

(4)  Assuming  that  the  initial  voltage  across  the  RC  network  up  is 
equal  to  zero,  due  to  the  zero  state  response  of  RC  circuit,  the 
dynamics  of  up  is 
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Fig.  2.  Battery  test  schedule. 


up(k+l)  =  e  At/Tup(k)  +  /battRp(l -e  At/T) 
r  =  KpCp  (7) 

up(0)  =  0 


where  k  is  the  time  step.  The  above  procedure  at  different  SOC 
values  is  repeated  to  establish  the  SOC  dependency  of  the  model 
parameters. 

3.2.  Battery  parameter  estimation  result 

The  modeling  results  include  the  following: 

(1)  The  Voc— SOC  curve  indicated  in  Fig.  5(a)  is  fitted  by  a  sixth- 
order  polynomial  with  95%  confidence  bounds  and  a  RMSE 
of  0.008328. 

(2)  The  relationship  between  Ro  and  SOC  in  Fig.  5(b)  is  fitted 
using  a  quadratic  polynomial  with  95%  confidence  bounds 
and  a  RMSE  of  0.0008608. 

(3)  As  shown  in  Fig.  5(c)  and  (d),  the  parameters  of  the  RC 
network  (Rp,  Cp)  as  a  function  of  SOC  are  approximated  by 
exponential  functions.  The  r  is  the  time  constant  of  RC 
network. 


Time/hour 


Time/hour 


Fig.  4.  Zoom-in  (120  s)  discharge  pulse  of  the  LiNMC  cell  at  round  90%  SOC  at  22  °C. 


3.3.  Parameters  estimation  result  discussion 

The  model  identified  was  built  in  Simulink  environment,  and 
the  evaluation  result  in  the  hybrid  pulse  test  is  shown  in  Fig.  6. 
There  is  good  agreement  between  the  model  and  the  experimental 
data  in  the  training  dataset  at  22  °C  (i.e„  hybrid  pulse  test),  and  the 
maximum  relative  error  is  less  than  1.5%. 

A  series  of  “unseen”  tests  were  conducted  to  verify  and  validate 
the  model.  The  robustness  of  the  model  is  validated  based  on 
different  dynamic  loading  profiles  at  22  °C,  as  shown  in  Figs.  7  and 
8.  The  mean  relative  errors  are  0.1617%  and  0.144%  in  the  DST  and 
FUDS  tests,  respectively.  The  maximum  relative  errors  in  both  tests 
are  less  than  1%. 

Fig.  9  shows  the  errors  of  the  nominal  model  at  different  cell 
aging  levels  under  22  °C.  It  indicates  that  as  the  aging  number  N 
increases,  the  mean  and  maximum  errors  of  the  nominal  model 
increase.  Based  on  the  quantitative  analysis,  the  model  parameters 
are  to  be  updated  when  a  significant  deviation  from  the  measured 
cell  voltage  is  detected.  Possible  trigger  variable  to  identify  a  loss  in 
model  accuracy  is  the  instantaneous  voltage  difference  from  the 
measured  value,  or  the  average  voltage  error  over  a  drive  cycle.  The 
actual  trigger  threshold  can  be  context  specific  and  relevant  to  the 
particular  criteria  of  the  target  application. 

Based  on  Fig.  9,  an  average  relative  model  error  of  0.7%  is  here 
chosen  as  the  trigger  threshold,  which  is  quite  useful  to  design  the 
time  scale  of  the  subsequent  SOH  estimator. 


Fig.  3.  Current,  voltage,  and  SOC  of  the  LiNMC  cell  in  the  training  dataset  at  22  °C. 
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(a)  Voc-SOC 


(b)RQ-SOC 
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(c)  R  -SOC 
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(d)  x-SOC 


O  Estimate  Result  - Fitness  Result  Pred  bnds  (95%) 


Fig.  5.  Parameter  estimation  result  in  HPPC  at  22  °C. 

4.  Combined  SOC  and  SOH  estimation 

The  increasing  estimation  error  results  from  the  variations  of 
model  parameters  as  the  battery  cell  ages.  To  keep  the  sufficient 


(a)  Terminal  Voltage  in  DST  at  22  °C 


3.2 

11 


11.48,  11.52  ,  11.56  , 

5  12  12.5  13 


13.5 


14.5 


(b)  Terminal  Voltage  Error  in  DST  at  22  °C 


Time/ Hour 


Fig.  7.  Validation  result  in  the  DST  test  at  22  °C. 


accuracy  of  SOC  estimation,  it  is  plausible  to  update  the  parameters 
through  the  combined  estimation  of  SOC  and  SOH,  whenever  the 
error  of  SOC  estimation  is  unacceptable.  Fig.  10  shows  the  multi¬ 
scale  state  estimation  framework.  By  offline  checking  the  accu¬ 
racy  decline  of  the  nominal  model,  we  can  determine  whether  it  is 
time  to  update  the  parameters  (i.e„  to  trigger  the  SOH  estimator). 
For  offline  parameters  recalibration  and  SOH  estimation,  a  fourth- 
order  EKF  is  used.  The  state  vector  is 


x  =  [SOC  Up  R0  l/CCap]T 

The  state  space  function  is  indicated  below, 


x(fc+l)=A(k)x(k)+B(k)/(k) 

y(k+l)=C(k)x(k) 

where 


0  0 


A(k) 
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Fig.  6.  Validation  result  in  the  hybrid  pulse  test  at  22  °C. 
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(a)Terminal  Voltage  in  FUDS  at  22°C 
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Fig.  8.  Validation  result  in  the  FUDS  test  at  22  °C. 
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Fig.  9.  Mean  and  maximum  relative  errors  of  the  nominal  model  at  different  battery  health  levels. 


Fig.  10.  Schematic  of  the  combined  SOC/SOH  estimation. 
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(a)  SOC  estimation  result 
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Fig.  11.  State  estimation  result  during  the  hybrid  pulse  test  in  the  cycle  22  at  22  °C. 


(a)  Estimated  capacity  result  at  N=22 


Fig.  12.  Capacity  estimation  with  different  EKF  initial  parameters  and  cycle  numbers  during  the  hybrid  pulse  test  at  22  °C. 
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Fig.  13.  Capacity  and  SOH  estimation. 


the  text  but  shown  in  Fig.  10.  As  introduced  in  Section  3,  the  average 
relative  error  threshold  r  is  0.7%. 

Fig.  11  shows  the  estimation  results  from  the  fourth  EKF  for  the 
fresh  battery  during  the  hybrid  pulse  test  in  the  cycle  22  at  22  °C. 
The  SOC  estimation  error  is  less  than  0.1%  during  the  most  of  the 
time,  and  the  estimated  resistance  is  closed  to  the  test  value.  The 
estimated  capacity  converges  toward  the  calibrated  value  reason¬ 
ably  whenever  the  initial  guess  starts  from  the  values  with  ±10% 
relative  error.  The  estimated  terminal  voltage  almost  superposes 
the  test  value,  and  the  mean  error  is  less  than  0.5%. 

Fig.  12  highlights  that  the  capacity  estimation  converges  toward 
the  test  value  when  the  battery  degrades  at  the  different  levels 
whenever  the  battery  is  almost  new  (N  =  22),  served  for  a  while 
(N  =  604)  or  used  heavily  (N  =  1273).  It  should  be  noted  that  the 
estimation  shows  good  convergence  when  the  covariance  param¬ 
eter  Po(4,4)  corresponding  to  the  capacity  estimation  is  0.5, 1,  5,  or 
10.  It  justifies  that  the  proposed  fourth  EKF  method  has  sufficient 
robustness  for  an  easy  implementation  in  practical  applications. 

4.1.  Results  discussion 


c(,i|  =  (e7  °) 


4.1.1.  Capacity  results  after  different  cycles 

The  battery  capacity  estimates  at  the  different  aging  levels  of  the 
battery  are  shown  in  Fig.  13.  The  estimates  with  red  (in  the  web 
version)  cross  markers  are  the  average  of  those  t  with  square  and 
triangle  markers.  It  is  clear  that  capacity  fade  is  predicted  with 
acceptable  accuracy.  SOH  defined  by  Eq.  (6),  is  also  shown  in  Fig.  13. 

SOH=  Ccap  x  100%  (9) 

^-Nominal 

where  Qap  is  the  actual  capacity  shown  in  Fig.  1  and  CNominai  is  the 
nominal  capacity  shown  in  Table  2. 


/(SOC)  =  a~i  ■  SOC6  +  a2  ■  SOC5  +  . .  ,a6  ■  SOC  +  a7. 

where  At  is  sample  time,  0.1  s.  Here, /(SOC)  is  the  fitted  polynomial 
of  OCV— SOC  relation.  EKF  algorithm  equations  are  not  elaborated  in 


4.1.2.  Resistance  results  after  different  cycles 

Fig.  14  shows  the  resistance  estimates  by  the  fourth  EKF  at  the 
different  battery  ages,  with  the  acceptable  accuracy  comparing  to 
the  ohmic  resistance  values  obtained  according  to  Eq.  (3). 
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Fig.  14.  Resistance  estimates  during  the  hybrid  pulse  test  in  the  different  cycles  at  22  °C. 
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Fig.  15.  Voltage  errors  before  and  after  model  recalibration. 

4.1.3.  Voltage  behavior  examination 

The  estimation  accuracy  for  the  battery  terminal  voltage  is  much 
improved  when  the  model  parameters  are  updated  by  means  of  the 
fourth-order  EKF,  especially  when  the  aging  of  the  battery  is  sig¬ 
nificant  (see  Fig.  15).  It  can  be  seen  that  the  voltage  errors  stay 
below  0.7%  after  updating  parameters. 

4.1.4.  SOC  estimation  results 

SOC  estimation  results  at  different  aging  levels  in  the  hybrid 
pulse  test  are  shown  in  Figs.  16—18.  Compared  to  the  second-order 
EKF  without  parameters  update,  the  offline  fourth-order  EKF  en¬ 
hances  the  accuracy  of  SOC  estimation  considerably  when  the 
battery  ages  significantly.  It  could  be  concluded  that  it  is  undesir¬ 
able  to  estimate  the  battery  SOC  without  considering  the  capacity 
and  resistance  recalibration  over  the  entire  battery  lifecycle.  Since 
the  fourth-order  joint  EKF  is  very  computationally  heavy  and  may 


be  subject  to  stability  issues,  the  SOH  estimator  in  this  paper  is 
intentionally  operated  offline,  and  the  related  capacity/resistance 
updates  are  instead  used  in  the  online  second-order  SOC  estimator, 
resulting  in  nearly  the  same  SOC  trajectory  as  that  in  the  offline 
fourth-order  one. 


5.  Further  discussion 

The  fourth-order  EKF  SOH  estimator  is  regularly  used  (e.g.,  after 
a  few  hundreds  of  cycles)  when  SOC  estimation  and  modeling  er¬ 
rors  become  unacceptable.  SOH  is  manifested  as  capacity  fading 
and  internal  resistance  increment  which  extremely  slowly  alter  in 
practice.  Thus,  unlike  the  SOC  EKF,  we  did  not  set  the  joint  fourth- 
order  EKF  online.  Offline  SOH  estimation  alleviates  the  requirement 
for  the  algorithm  stability  (i.e.,  we  can  use  sufficient  batch  of  data  to 
slowly  converge  to  accurate  SOH  estimates,  see  Figs.  11  and  12). 
Furthermore,  as  the  classic  book  [39]  indicates,  this  joint  fourth- 
order  EKF  is  an  adaptive  structure,  which  may  be  subject  to  sta¬ 
bility  problem.  Theoretically,  it  is  still  very  difficult  to  derive  a  proof 
of  stability  for  this  type  of  estimator.  However,  Ref.  [39  also 
demonstrates  that  if  model  parameters  are  slowly  varying 
adequately,  joint  EKF  is  stable  in  most  cases.  In  our  estimation 
scenario,  the  SOH  parameters  (i.e.,  capacity  and  internal  resistance) 
fulfill  this  slowly-varying  requirement  very  well,  which  has  a 
marginal  adverse  influence  on  the  EKF  stability,  leading  to  a 
convergent  estimator.  Please  notice  that  the  joint  estimator  is  off¬ 
line,  instead  of  real-time  realization  such  that  large  amounts  of  data 
can  be  exploited  to  assure  its  asymptotical  convergence.  After 
updating  the  batteiy  SOH  (the  model  parameters)  offline,  the  on¬ 
line  SOC  EKF  are  provided  with  constant  model  parameters  upda¬ 
ted,  and  thus  the  single  SOC  EKF  is  decoupled  with  the  fourth-order 
SOH  EKF.  This  is  essentially  different  from  two  online  joint  EKFs. 
Given  constant  parameters,  the  SOC  EKF  often  can  practically  meet 
stability  and  convergence  requirements. 


(a)  SOC  Estimate  Result 


Fig.  16.  SOC  estimation  result  for  the  fresh  battery  in  hybrid  pulse  test  (N  =  22,  22  °C). 
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Fig.  17.  SOC  estimation  result  for  the  moderately  aged  battery  in  hybrid  pulse  test  (N  =  604,  22  °C). 


(a)  SOC  Estimate  Result 


(c)  SOC  estimation  improvement  by  4th  order  EKF 


Fig.  18.  SOC  estimation  result  for  the  heavily  aged  battery  in  hybrid  pulse  test  (N  =  990,  22  °C). 


6.  Conclusions 

Recursive  least  squares  algorithm  is  first  employed  to  identify 
parameters  of  the  first-order  RC  model  for  a  LiNMC  battery  over 
its  cycle  lifespan  in  this  article.  Then,  the  performance  decay  of 


the  nominal  model  is  quantified  for  different  battery  aging  levels. 
Based  on  the  modeling  results,  a  multi-scale  observer  approach 
is  proposed  to  implement  combined  SOC/SOH  estimation  over 
the  lifespan  of  the  lithium-ion  battery.  The  SOC  is  estimated  in 
real-time  by  using  a  second-order  EKF,  and  the  SOH  (the  capacity 
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and  internal  ohmic  resistance)  is  updated  offline  in  a  fourth- 
order  EKF.  The  trigger  time  of  the  SOH  estimator  is  carefully 
specified  based  on  the  quantified  model  accuracy  deterioration 
with  increasing  battery  aging.  Extensive  experimental  data  are 
applied  to  demonstrate  that  the  developed  battery  monitoring 
algorithm  can  provide  very  accurate  SOC  and  SOH  estimates 
without  heavy  computational  burden  and  occurrence  of  insta¬ 
bility/divergence. 
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